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I. INTRODUCTION 

Simple discrete models are valuable tools for exploring 
phase transitions, both in equilibrium and out of equi- 
librium. A paradigmatic example of the first kind is the 
kinetic Ising model (heat-bath dynamics), which can be 
formulated as a Probabilistic Cellular Automaton (PCA) 
(see Appendix and Ref. [1]). One of the simplest ex- 
amples of a system exhibiting out of equilibrium phase 
transitions is the Domany-Kinzel (DK) model [2], which 
is a natural extension of the directed percolation (DP) 
problem in the language of cellular automata. 

Frustrations play a central role in characterizing the 
phase diagrams of many simple models. One can have 
geometric frustrations, like in a triangular antiferromag- 
netic Ising model, or fluctuating coupling, like in spin 
glasses [3] or p-spin [4] models. For instance, in the p-spin 
model, the character (ferro or anti-ferro) of the interac- 
tion varies widely with the local magnetization (i.e. even 
a flip of a single spin may invert the interaction). We 
discuss here a simpler one dimensional model in which 
frustrations play a central role. This model is a variation 
of the droplet model in which the particles tend to repel 
each other when they are dispersed but in which clus- 
ters, once formed, cannot break and only particles near 
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the surface can eventually leave them. 

The model originated from the schematization of the 
mechanism of opinion formation in a homogeneous (no 
leader) society [5], in which it was assumed that a co- 
herent local community exerts a very strong social pres- 
sure on an individual's opinion. However, there is the 
possibility of disagreement with a weak local majority 
depending on the individual's education. In the case of 
a diffuse non-conformistic attitude, people tend to act in 
the opposite way of the local majority (antiferromagnetic 
coupling), introducing frustrations. The presence of ab- 
sorbing states may cause the formation of large coherent 
clusters, whose interactions give way to interesting pat- 
terns. 

This model could also be applied to a system of 
charged, magnetized metallic particles. These repel each 
other due to their charge, but are attracted due to their 
magnetic coupling. With a proper choice of parame- 
ters, one may have repulsion when particles are scattered. 
However, once a cluster of particles in contact appears, 
the mobility of charges may lower the electric field among 
this cluster and an external particle, and the interaction 
may become attractive. Other examples can be found in 
aggregation models. 

In what follows, we discuss a one dimensional spin 
model characterized by two coupling factors; one that be- 
haves like the Ising (magnetic) coupling, i.e., for a given 
spin, grows linearly with the local field, and a nonlinear 
term that is very low for a weak local field, but grows 
quickly when the local field exceeds a certain threshold. 

In one dimension, no true phase transition appears for 
finite-range and finite-strength coupling. On the other 
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hand, models presenting absorbing states, like the DK 
one, do exhibit noncquilibrium phase transitions. The 
presence of absorbing states may be related to the diver- 
gence of some coupling, that becomes infinite (see Ap- 
pendix and Rcf. [1]). Thus, in our model, we simply as- 
sume a two-level nonlinear coupling: zero for local field 
below a given threshold (corresponding to the parameter 
Q in the following) and infinite otherwise. Thus, in to- 
tal we have four parameters: the range R, the external 
magnetic field H and local coupling J like in the Ising 
model, and the threshold Q. We denote these models by 
RQ, indicating the range and the threshold, say i?3Ql, 
R5Q2, ... 

The model can also be discussed as a one dimensional 
probabilistic cellular automaton. In this guise, it may 
be considered an extension of the DK model, or, more 
precisely, of the R3Q1 model [6]. This model exhibits a 
richer phase diagram than the DK one, and we shall show 
that for R — > oo, the RQ model presents novel features, 
like a disorder-chaotic transition. 

It is not surprising that frustrations gives rise to dis- 
ordered behavior, which may be sometimes considered 
chaotic from a microscopic point of view. We investigate 
this aspect using an original chaotic observable, corre- 
sponding to a finite-distance Lyapunov exponent. How- 
ever, when renormalized (or locally averaged), a "disor- 
dered" behavior simply corresponds to a stable almost- 
zero magnetization. In the mean-field approximation, 
it corresponds to a fixed point of observables. In our 
model also a coherently-chaotic phase appears, in which 
large patches of sites oscillate widely, corresponding to a 
chaotic map in the mean-field approximation. 

In summary, we observe that the usual transitions from 
ordered into active phases becomes much richer, being 
preceded by a transition from coherently chaotic to sim- 
ply chaotic, then to more irregular states, the appearance 
of a parameter-insensitive disordered phase and finally an 
"escape" to ordered (quiescent) states with a discontinu- 
ous (first-order) character. The origin of this behavior is 
the competition between "linear" antiferromagnetic and 
"nonlinear" ferromagnetic couplings. 

The paper is organized as follows. In Sec. II we present 
the RQ model both as a spin system and a PCA. In 
Sec. Ill we present numerical results of the model begin- 
ning with a review of the simplest non trivial case R = 3, 
Q = 1 and Sec. IV is devoted to a mean field discussion 
of the model. The paper ends with some conclusions. 



II. THE MODEL 

Both the kinetic Ising model (see Appendix) and the 
DK one are Boolean, one dimensional, R = 2 PCA in 
which the state of a given cell depends probabilistically 
on the sum S of the states of its two nearest neighbors 
at the previous time step. That is, there are three dif- 
ferent conditional probabilities, t(1\S), S = 0,1,2 that 
the future state of the cell is 1 given S neighbors in state 



1. In the DK model r(l|0) = which means that there 
is an absorbing phase corresponding to the configuration 
where the state on every cell is zero. The presence of 
more than one absorbing state can induce different be- 
haviors and trigger the appearance of universality classes 
different from the usual directed percolation (DP) one [7- 
9]. 

We denote by s\ G {0, 1} the state at site i = 0, . . . , L — 
1 and time t = 0, 1 . . . . All operations on spatial indices 
are assumed to be modulo L (periodic boundary condi- 
tions). The range of interactions is denoted by R, that 
is, s* +1 depends on the states at time t in a neighbor- 
hood Mi that contains the R nearest neighbors of site 
i. It is convenient to introduce also the spin notation: 
a\ = 2s\ - 1, g\ e {-1, 1}. 

The spin a\ "feels" the local field V{m\) = H+[J + 
K{m\)]m\, where H is an external field, J a coupling 
constant and the local magnetization is defined as 



E-5 



2S\ - R, 



where 
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The presence of absorbing states is due to the nonlinear 
coupling K given by 



K{m) = K(m(S,R),Q) 



if S < Q, 

if S > R - Q, 

otherwise. 



with k a constant. This term is relevant only if it is 
dominant with respect to the linear one. We thus choose 
the limit k — > oo, so that the transition probabilities 
t(1\S) are given by 



r(l\S) 



1 + cxp {-2 [H + Jm{S, R)]} 



if S < Q, 
if S > R - Q, 

otherwise. 

(1) 

For Q = we recover the usual heat-bath dynamics 
for an Ising-like model with reduced Hamiltonian Tt = 
J2iV( m i)&i- For R = 3, Q = 1 we have the model of 
Ref. [6] with two absorbing states. The quantities H 
and J range from — oo to oo. For easy plotting, we use 
j = [1 + cxp(-2J)]" 1 and h = [1 + cxp(-2 J ff)]" 1 as 
control parameters, mapping the real axis (— oo, oo) to 
the interval [0,1]. 

The fraction c of ones in a configuration and the con- 
centration of clusters p are defined by 
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Both the uniform zero-state and one-state correspond to 
p — > in the thermodynamic limit, while the active state 
corresponds to p > 0. 
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FIG. 1: Phase diagram of R3Q1 model. The upper-left region 
corresponds to the absorbing state 0, the upper-right region 
to the absorbing state 1 and the lower region to the disordered 
phase. In the inset we show the asymptotic value of p as a 
function of j for H = (h = 0.5). 




FIG. 2: Typical space time patterns of the R3Q1 model for 
H = (h = 0.5) with L = 256. (A - top left) Chaotic de- 
terministic pattern j = 0.03125. (B - top right) Disordered 
phase j — 0.5. (C - bottom left) competition between quies- 
cent phases j — 0.54688 (first-order transition). 



III. NUMERICAL RESULTS 

Let us first review the i?3Ql case, whose phase dia- 
gram is reported in Fig. 1. A first order phase transition 
occurs on the vertical line that ends at a bicritical point 
where it meets two second order phase transition curves. 
In the neighborhood of the first order phase transition 
line a typical pattern of the system is composed by large 



patches of zeroes and ones, separated by disordered zones 
(walls) whose width does not grow in time. These walls 
perform a sort of random motion and annihilate in pairs 
(see Fig. 2-C). By lowering J an active, disordered phase 
appears (see Fig. 2-B). The transition between the or- 
dered (absorbing) and disordered (active) phases occurs 
by destabilization of the width of the walls, that perco- 
late in the whole system. In the limit J = — oo (j = 0) 
we have typical "chaotic" class-3 cellular automata pat- 
terns [10], as shown in Fig. 2- A. Damage spreading analy- 
sis [6] show that inside the active phase there is a region in 
which a variation in the initial configuration can influence 
the asymptotic configuration. This region may be called 
chaotic. However, in the mean field approximation, this 
region simply corresponds to a stable fixed point of the 
average density. Indeed, the fluctuations of the density 
remains quite small and a spatial coarse-graining would 
generate a simple random pattern. 

This underlining "disorder" is revealed also by the 
asymptotic value of p for H = (h = 0.5) as a function 
of j, shown in the inset of Fig. 1. It exhibits a monotonic 
behavior, with a maximum for j = (corresponding to 
the smallest average length of homogeneous clusters) and 
a continuous phase transition at the bicritical point, that 
separates the active (p > 0) and quiescent (p — 0) phases. 

By increasing R, new features appear. Typical pat- 
terns are shown in Fig. 3 for different values of j. As 
illustrated in Fig. 4, p is no longer a monotonic function 
of j, and a new, less disordered region appears inside the 
active one for small values of j. The transition between 
the active and the quiescent phases become sharper with 
increasing R, as shown in Fig. 4. 

The pattern shown in Fig. 3-a is reminiscent of 
"chaotic" deterministic cellular automata. We refer to 
this region as coherently- chaotic, since it corresponds 
to "irregular" coherent oscillations of large patches of 
sites. This region of the parameters is dominated by an 
almost-deterministic behavior, and the presence of many 
metastable "absorbing" states, revealed by transient reg- 
ular patterns in Fig. 3-a. As we discuss in Sec. IV, the 
mean field analysis gives a chaotic map for large anti- 
ferromagnetic values of J, but we were unable to find 
a clear order parameter that numerically distinguishes 
this region from the broader "chaotic" one in which the 
asymptotic configuration exhibits a dependence on varia- 
tions in the initial configuration. In Sec. Ill A, we analyze 
the chaotic region by means of (finite size) Lyapunov ex- 
ponents. 

By increasing J, the coherent patches shrink and the 
Lyapunov exponent decreases and finally becomes nega- 
tive. The system asymptotically loses its dependence on 
the initial conditions and is dominated by the stochastic 
components. This phase is termed irregular, and appears 
completely random. Fig. 3-b shows a typical pattern 
at the boundary between the chaotic and the irregular 
phase. 

Our simulations were carried out using the fragment 
method [11], in which a set of configurations (replicas) 
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FIG. 3: Typical space time patterns (L = 256) for model 
R11Q1 and H = 0, starting from a disordered initial condi- 
tion with p = 0.5. (A) [top left, j = 0.056250]: "coherently- 
chaotic" phase. One can observe rare "triangles" with "base" 
larger that R, corresponding to the unstable absorbing states, 
and other metastable patterns, corresponding to additional 
absorbing states for J —* — oo. (B) [top right, j = 0.421875]: 
At the boundary between the chaotic and the irregular phase, 
the only locally absorbing states are (rare) patches of and 
1. Most of the pattern looks random, with some "triangles" 
reminiscent of chaotic CA. (C) [bottom left, j = 0.478125]: 
Disordered phase, the pattern looks random and the differ- 
ence between patterns obtained with different values of the 
parameters (and the same random numbers) is vanishing. (D) 
[bottom right, j — 0.562500]: Quiescent phase. In this phase 
the only stable states are the absorbing ones. The boundaries 
separating the phases move randomly until coalescence. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

j 

FIG. 4: Behavior of p for H = 0. RSQ1 (long dashed), M1Q1 
(solid), i?41Ql (dashed) and i?81Ql (dotted). Letters corre- 
spond to patterns in Fig 3. 



are updated with different values of parameters using the 
same random numbers, and the same initial configura- 
tion. 

By observing in rapid sequence the patterns generated 
with increasing values of J, one observes a sudden "freez- 
ing" of the (random) image, just before the transition 
between the active to quiescent phases, corresponding to 
Fig. 3-c. This effect is due to the insensitivity of the 
patterns not only to the initial condition, but also to the 
parameters, see Sec IIIB. We term this phase disordered. 

Finally, the quiescent phase is asymptotically domi- 
nated by the two absorbing states, with the usual anni- 
hilating walls dynamics, Fig. 3-d. 

Fig. 4 shows that the active-quiescent transition be- 
comes sharper when R is increased, and that the clus- 
ter density exhibits a cusp at the transition. By en- 
larging this region one sees that for H = (Fig. 5) 
the cluster density exhibits two sharp bends, while for 
H = 0.42, (Fig. 6), only one bend is present. This is 
reminiscent of the universality class change (parity con- 
servation, DP) for H = 0, H y£ 0, respectively, in the 
R3Q1 model [6]. Notice also that the irregular-disorder 
transition at j ~ 0.494 is rather peculiar, since p first 
decreases and then suddenly jumps to high values. 

The bends are not finite size or time effects. As shown 
in Fig. 3-c, in this range of parameters the probability of 
observing a local absorbing configuration (i.e. patches of 
zeroes or ones) is vanishing. All other local configurations 
have finite probability of originating zeroes or ones in the 
next time step. The observed transitions are essentially 
equivalent to those of an equilibrium system, that in one 
dimension and for short-range interactions cannot exhibit 
a true phase transition. The bends thus become real 
salient points only in the limit R — > oo. 



Chaos and Lyapunov exponent in cellular 
automata 



State variables in cellular automata are discrete, and 
thus the usual linear stability analysis classifies them as 
stable systems. The occurrence of disordered patterns 
and their propagation in stable dynamical systems can 
be classified into two main groups: transient chaos and 
stable chaos. 

Transient chaos is an irregular behavior of finite life- 
time characterized by the coexistence in the phase space 
of stable attractors and chaotic non attracting sets, 
namely chaotic saddles or repellers [12]. After a tran- 
sient irregular behavior, the system generally collapses 
abruptly onto a non-chaotic attractor. This kind of be- 
havior is reminiscent of some class-4 (nonchaotic) de- 
terministic cellular automata (DC A) [10], and can be 
present in continuous systems having a discrete dynamics 
as a limiting case [13]. 

Stable chaos constitutes a different kind of transient 
irregular behavior [14, 15] which cannot be ascribed to 
the presence of chaotic saddles and therefore to diver- 
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FIG. 5: Comparisons between numerical (thin line) and mean- 
field (thick, dotted line) results for R81Q1 and H = 
(h = 0.5). The estimated critical values are j* ~ 0.493827 
and jii ~ 0.51384. This figure has been obtained with much 
larger simulations of the corresponding line in Fig. 4, and 
more details emerge. The patterns shown in Fig. 3 (which 
however refer to R — 11) illustrate the typical behavior of the 
system to the left (B), in the middle (C) and to the right (D) 
of the plateau. 




FIG. 6: Comparisons between numerical (thin line) and mean- 
field (thick, dotted line) results for i?201Q5 and H = 0.42 
(h = 0.7). The two set of pluses mark the period-two re- 
gion. The estimated critical values are j* ~ 0.49740 and 
ju ~ 0.50369. 



gence of nearby trajectories. Moreover, the lifetime of 
this transient regime may scale exponentially with the 
system size (supertransients as defined in Refs. [14, 15]), 
and the final stable attractor is practically never reached 
for large enough systems. One is thus allowed to assume 
that such transients may be of substantial experimental 
interest and become the only physically relevant states 
in the thermodynamic limit. 

The emergence of this "chaoticity" in class-3 (chaotic) 
DC A dynamics, is effectively illustrated by the damage 
spreading analysis [16, 17], which measures the sensitiv- 



FIG. 7: Contour plot of the maximum Lyapunov A exponents 
for different values of R and j, for H = 0. The solid line 
represent the boundary between the A > phase and the 
A = — oo one. 



ity to initial conditions and for this reason is considered 
as the natural extension of the Lyapunov technique to 
discrete systems. In this method, indeed, one monitors 
the behavior of the distance between two replicas of the 
system evolving from slightly different initial conditions, 
or the speed of propagation of a disturbance [18]. The 
dynamics is considered unstable and the DCA is called 
chaotic, whenever a small initial difference between repli- 
cas spreads through the whole system. On the contrary, 
if the initial difference eventually freezes or disappears, 
the DCA is considered non chaotic. 

Due to the limited number of states of the automa- 
ton, damage spreading does not account for the maximal 
"production of uncertainty", since the two replicas may 
synchronize locally just by chance (self-annihilation of 
the damage). Moreover, there are different definitions of 
damage spreading for the same rule [19]. 

To better understand the nature of the active phase, 
and up to what extent it can be denoted chaotic, we 
extend the finite-distance Lyapunov exponent defini- 
tion [20] to probabilistic cellular automata. A similar 
approach has been used in Ref. [21], calculating the Lya- 
punov exponents of a Kauffman random Boolean net- 
work in the annealed approximation. As shown in this 
latter paper, this computation gives the value of the (clas- 
sic) Lyapunov exponent obtained by the analysis of time- 
series data using the Wolf algorithm [22]. 

Given a Boolean function f(x,y,...), we define the 
Boolean derivative df/dx with respect to x by 

5/ fl iff(\x-l\,y,...)^f(x,y,...), 
dx [0 otherwise, 

which represents a measure of sensitivity of a function 
with respect to x. The evolution rule of a probabilistic 
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FIG. 8: Susceptibility % (continuous line) and average number 
of clusters p (dashed line) for R11Q1, H = 0, Aj = 0.01 



cellular automaton can be written as 



r < r(l\Sl), 
otherwise 



where r is a random number uniformly distributed in 
[0, 1]. The Boolean derivative with respect to Xj is eval- 
uated by using the same r in the comparison implied in 
Eq. (2). 

For a PCA, we can thus build the Jacobian matrix 
Jij = dxl +1 /dxj, i,j = 0, ...,L — 1. The maximum 
Lyapunov exponent A can now be defined in the usual 
way as the expansion rate of a tangent vector v(t), whose 
time evolution is given by 

v(t + 1) = Jv(t). 

When all the components of v become zero, A = — oo 
and no information about the initial configuration may 
"percolate" as t — > oo [20]. This maximum Lyapunov 
exponent is also related to the synchronization properties 
of cellular automata [23]. 

A preliminary numerical computation of A for our 
model is reported in Fig. 7. It can be noticed that the 
boundary j c (R) of A > is not monotonic with _R, reach- 
ing a maximum value for R ~ 11. By comparisons with 
Fig. 4, one can see that the chaotic phase is included in 
the irregular one. 



B. Parameter dependence 

We can numerically characterize the irregular-disorder 
transition by looking at the sensitivity of space-time pat- 
terns with respect to the variations of parameters. Using 
the fragment method [11] a set of configurations (repli- 
cas) are updated with different values of parameters using 
the same random numbers, i.e. with the same disordered 
field and initial conditions. The quantity 



x(J) 



= lim lim } 



4(J + AJ) - 8*(J)\ 



is a measure of the susceptibility. For large correlations 
one expects very small differences among replicas. As 
shown in Fig. 8, this susceptibility is strongly correlated 
to the average number of clusters p. 



IV. MEAN FIELD APPROXIMATION 

In the simplest mean field approximation, the density 
c evolves in time according to 



c(f + 1) = W^ s) (c(t)) T {l\S), 



(3) 



with 



w^ s \c)=Qy s (i- C ) R - s . 

An example of a mean-field phase diagram is reported in 
Fig. 9. 

For large R, W^ R ' S ^ can be approximated by 



y/2nc{l - c)R 



exp 



R(S/R - cf 
2c(l - c) 



(4) 



and the summation can be replaced by an integral. The 
parameters of the resulting equation may be rescaled by 
using J = JR and Q = Q/R. In the limit R — > oo (J 
and Q fixed), c(t + 1) = f(c(t); H, J) with 



1 + cxp(-2(# + J(2c - 1)) 



if c < Q, 
if c > 1 - Q, 

otherwise. 



(5) 



The requirement of constant rescaled variables means 
that the variations of J for large R triggers large vari- 
ations of J, and consequently the transition region be- 
comes sharper. 

This approximation is supposed to be particularly 
good in the disordered phase, i.e. for J, H ~ 0, as shown 
in Fig. 5. In the finite-dimensional case, one has to add 
a noise term, representing the influence of neighboring 
sites, whose amplitude is of the order of the width of the 
Gaussian in Eq. (4), \/ c(l — c)/R. 

In this mean-field approximation, the order parameters 
c and the numerical density of clusters p are related by 
p = 2c(l - c). 

For J, H ~ the mean field approximation gives a 
stable fixed point c*(H, J) ^ 0, 1 for the density, which 
assumes the maximum value c* = 0.5 for J = H = 0. 

The origin of the irregular-disordered and disorder- 
quiescent phase transitions marked by (I) and (II) re- 
spectively in Figs. 5 and 6 is due to the loss of stability 
of the mean-field fixed point c* ^ 0, 1 (see Eq. (5) and 
Fig. 10). 
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FIG. 9: Mean-field phase diagram for R — 81. The upper 
curves correspond to the transition II from the quiescent to 
the disordered phases, R81Q5 (dashed) and R81Q1 (dotted). 
The solid curve corresponds to the transition I from the dis- 
ordered to the irregular phase. 
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FIG. 10: Mean field map Eq. (3) in the neighborhood of tran- 
sition I (J ~ —1, dashed line) and II (J ~ 1, solid line), for 
H = and Q = 0.02 



The transition I is given by 



ct = f(4;Hj); 



df(c;H, J) 



dc 



(6) 



which would induce a period-doubling cascade in the 
mean-field approximation (disregarding the presence of 
absorbing states). The solution of Eq. (6) is 

4 = ~ (i + Vi + J- 1 ) , 



9 l0g I -T 



The critical value Jj" can be obtained numerically once 
given the value of H* = H*(J*), and corresponds to the 
left bend in Fig. 5 and 6. 



For H = (Fig. 5) the period-doubling instability 
brings the local configuration into an absorbing state, 
and the lattice dynamics is therefore driven by interac- 
tions among patches which are locally absorbing. This 
essentially corresponds to the dynamics of a Determin- 
istic Cellular Automata (DCA) of chaotic type, i.e. a 
system which is insensitive to infinitesimal perturbations 
but reacts in an unpredictable way to finite perturbations 
[13]. In other words, in this region the original stochastic 
model behaves like a "chaotic" deterministic one after a 
coarse- graining of patches. We expect that this corre- 
spondence will become more and more exact with grow- 
ing R. From a theoretic-field point of view this means 
that the renormalization flux tends towards a "chaotic" 
model instead of the usual fixed-point dynamics. 

For H = 0.42, as shown in Fig. 6, the period-two phase 
has a finite amplitude, before falling into the DCA-like 
dynamics by reducing J. 

For J > (transition II) we have 



f{ct l ;H,J) = 



and thus the critical value Jj*j is 




J* 
■hi 



1 



2Q- 1 



±H 



log 



1-0 



Once that a portion of the lattice has been attracted to 
an absorbing state, it pulls the neighboring regions to this 
same state, due to the ferromagnetic (J > 0) coupling. 

This approximation, disregarding fluctuations, overes- 
timates the critical value Jfj, as shown in Figs. 5 and 
6. 

We analyze extensively the mean-field behavior of the 
system for H = and J = — oo. As shown in Fig. 11, for 
a given value of R, there is always a critical Q c value of 
Q for which the active phase disappears, with an approx- 
imate correspondence Q c ~ 2/5R. One can also observe 
that the chaotic oscillations of the mean-field map, that 
roughly corresponds to the coherent-chaotic phase of the 
model, appear only for large values of R. Since the ab- 
sorbing states are always present, there chaotic oscilla- 
tions may bring the system into the quiescent phase, as 
shown by the "hole" to the right of triangles in Fig. 11. 

The correspondence between the chaotic behavior of 
the mean-field approximation and the actual behavior of 
the system will be the subject of a future work. 



V. CONCLUSIONS 

We have investigated a general one dimensional model 
with extended-range interactions and symmetric absorb- 
ing states. The model is characterized by a competing 
ferromagnetic linear coupling and an antiferromagnetic 
nonlinear one. 
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FIG. 11: Case H = 0, J = — oo: the mean-field R — Q phase 
diagram of the mean-field approximation. Points mark pa- 
rameter values for which the absorbing states are the only 
stable attractors. A plus sign denotes period-2 temporal oscil- 
lations between absorbing states, a star denotes the presence 
of a stable point at c = 0.5, a cross (circle) denotes period-two 
(four) oscillations between two non-zero and non-one densi- 
ties, triangles denote chaotic oscillations. 



By means of numerical simulations and mean field 
approximations we have shown that a chaotic phase 
is present for strong antifcrromagnetic coupling. This 
phase may be identified as a stable-chaotic region, in 
which the behavior of the (originally stochastic) system is 
essentially deterministic and its behavior is highly irreg- 
ular and essentially unpredictable. The mean field map 
exhibit a chaotic behavior for a large interaction range R, 
and this behavior is reflected in the appearance of many 
metastable states in the system, for extremely strong an- 
tifcrromagnetic coupling. 

A disordered phase, insensitive of parameter varia- 
tions, appears at the boundary between the active and 
the quiescent ones, and the transitions appear to be of 
equilibrium type, i.e. truly salient points only in the limit 
R — > oo. 

We expect that in higher dimensions one can recover 
some aspects of these phase transitions without impos- 
ing the presence of absorbing states, i.e. using finite cou- 
plings. 
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APPENDIX A: EQUIVALENCE BETWEEN 
DYNAMIC ISING MODEL AND CELLULAR 
AUTOMATA 

As an illustration, we present here a derivation of the 
equivalence of the kinetic Ising model (here we choose 
heat bath) with a cellular automaton (the Domany- 
Kinzel model) in order to elucidate the role of infinite 
coupling parameters and absorbing states. This deriva- 
tion is similar to that of Ref. [1] but more general. 

The Ising model is defined by the couplings among 
spins. The configuration at time t is denoted as er = 
o"i , <T2 , . . . and the configuration at time t + 1 as er' = 

a[ , <j' 2 , Let us write the temperature-dependent 

Hamiltonian as 

^( <T ) = ^ #(ot-i)OtO»+i)) 

i 

with 

H(x, w, y) = J {a) w + J w xw + J {2) wy + J {3) xwy. (Al) 

The transition probabilities r must obey the detailed 
balance condition 



r(cr\*>) 



exp(H(cr') - H(a)). 



In each step we can update in parallel all even or odd 
sites, obtaining 

T ( tT '\ <T ) _ TT T~((Ti-i, gj, gj + i|g-j_i,(Tj,(Tj + i) 
T((T\(T') A.A r(cr i _l,(Tj,C7 i+ l|c7i_l,Cr-,(Tj + i) ' 

where the product is restricted to either even or odd sites. 
The detailed balance condition can thus be satisfied lo- 
cally. Choosing the Heat Bath dynamics 



t(x, l,y\x,w,y) 



exp(-H(x, 1,2/)) 



exp(-H(x, -1, y)) + exp(-H(x, 1, y)) 



the transition probabilities does not depend on the 
present value of the spin w — a\ so that the lattice 
(with even or infinite lattice sites) may be decoupled into 
two noninteracting sublattices with the same geometry of 
the DK model. From now on, to be coherent with the 
usual cellular automaton notation, we shall express the 
transition probabilities in terms of the Boolean variables 
Si = {<7i + l)/2, and we shall denote the local field as 
h(a,b) = H(2a- 1,1,26- 1). 

Let us denote the DK transition probabilities 
r(s-|si_i,s i+ i) as 



t(1|00) = £, t(1|01) = t(1|10) = p, 
We get for the heat bath dynamics 

1 



r(l|ll) = (jr. 



r(l|o6) = 



1 



exp( 



-2h{a,b)) 
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and thus Finally, we have 

u/ u\ 1 ^ r(l\ab) (0) 1 e 1 q 1 p 
h(a, b) = - In --— . = - In + - In —2— + - In — l —, 

2 l-r(l|a6) 8 1-e 8 l-q 4 1-p 

Substituting into Eq. (Al) one obtains a linear system = = — - In - — - + - In ^ - 

J(0 ) _ J( i) _ J( 2) J(3 ) = I ln _f_ j(3) = I In + 1 In 1 In 

2 1 — £ 8 1 — e 8 1— 4 1— p 

j(0) + j(i) _ j(2) _ j(3) = 1 i n _P_ i n the limit e -> all couplings become infinite. 

2 1 — p 

J(0)_ J(l) + J(2)_ J(3) = l ln . P 



2 1 -p 

J(0) +J (l) + J (2) + J (3) = l ln . ? 



2 l- g 
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